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Several heavy ion experiments at SPS have measured azimuthal distributions of particles with 
respect to the reaction plane. These distributions are deduced from two-particle azimuthal correla- 
tions under the assumption that they result solely from correlations with the reaction plane. In this 
paper, we investigate other sources of azimuthal correlations: transverse momentum conservation, 
' ' which produces back-to-back correlations, resonance decays such as A ^ pvr, HBT correlations and 

I final state interactions. These correlations increase with impact parameter: most of them vary with 

■ the multiplicity A'^ like 1/A'^. When they are taken into account, the experimental results of the 

' NA49 collaboration at SPS are significantly modified. These correlations might also explain an 

important fraction of the pion directed fiow observed by WA98. Data should be reanalyzed taking 
into account carefully these non-fiow correlations. 
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I. INTRODUCTION 



' In ultrarelativistic heavy ion collisions, azimuthal distributions of the produced particles with respect to the reaction 
^ ,. plane are of particular interest. Their anisotropics, which we refer to as "flow", result from the interactions between 
' the produced particles. They thus probe the hot stages of the collision and may provide valuable information on the 
[ state of the interaction region: thermalized or not, equation of state, etc... In the recent years, flow has been 

■ measured by many experiments at the Brookhaven AGS with Au projectiles [|[-|| and at the CERN SPS with Pb 
' projectiles p|-p^. Since the orientation of the reaction plane is not known a priori, flow measurements are usually 

■ extracted from two-particle azimuthal correlations. This is based on the idea that azimuthal correlations between two 
particles are generated by the correlation of the azimuth of each particle with the reaction plane. 

I The assumption that this is the only source of two particle azimuthal correlations dates back to the early days of 
"'Y' ' the flow It still underlies the analyses done at ultrarelativistic energies, both at AGS and at SPS (see, however, 
Ip^). However, there are various sources of direct azimuthal correlations between particles, which do not involve the 
reaction plane. We have shown in a recent paper [|l^ that the quantum correlations due to the HBT effect yield 
^ ' azimuthal correlations. When they are taken into account, the results obtained for the charged pion flow by the NA49 
collaboration ^ change significantly. Here, we investigate other sources of non-flow azimuthal correlations. A first 
source is the condition that the total transverse momentum of the outgoing particles is zero: this gives a back-to-back 
correlation between their momenta, which can be evaluated quantitatively. This effect is usually taken into account 
^ ' at lower energies but has been neglected so far at AGS and SPS. Other correlations are due to resonance decays 
I such as A ^ piT and p — > tttt. Finally, Coulomb and strong interactions between pairs of particles with low relative 
velocities (which we refer to as "final state interactions" ) produce small angle azimuthal correlations . 

Most of these correlations are typically of order 1/N, where N is the number of particles emitted in the collision. 
At SPS energies, correlations due to fiow are so small that additional correlations become of the same order, and 
cannot be neglected. The 1/iV dependence also determines the variation of non-fiow correlations with the centrality 
of the collision. Therefore they increase with impact parameter up to very peripheral collisions; on the other hand, 
they do not vanish for central collisions where the fiow is zero by symmetry. They are very important at large impact 
parameters and should be taken into account, in particular, when studying the centrality dependence of the flow, 
which has been recently proposed as a sensitive probe of the phase transition to the quark-gluon plasma pO|-p^ . 

In section ||, we briefly recall how the flow can be extracted from two-particle azimuthal correlations. We use the 
same method and notations as in Ifl7| . The various sources of non-flow correlations are reviewed and their effects are 



estimated in section [II. The experimental data obtained by the NA49 and WA98 collaborations for Pb-Pb collisions 
at SPS are discussed in section We show in particular that the effect of momentum conservation alone is large 
enough to reverse the sign of the proton directed flow measured by NA49 which is used to define the reaction 
plane. A decays produce correlations which are of the opposite sign, and of the same order of magnitude, although 
we do not attempt to calculate them accurately. For WA98, final state interactions play an important part, and may 
explain a large fraction of the observed pion directed flow. Our conclusions are given in section 
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II. TWO-PARTICLE AZIMUTHAL CORRELATIONS AND STANDARD FLOW ANALYSIS 



In this section, we define the flow and non-flow contributions to two-particle azimuthal correlations, and we show 
how azimuthal distributions with respect to the reaction plane can be extracted from these correlations. 

We call "flow" the azimuthal correlations between the outgoing particles and the reaction plane. These are best 
conveniently characterized in terms of the Fourier coefhcients f„ psf which we now define. We choose a coordinate 
system in which the x axis is the impact direction, and (a;, z) the reaction plane, cf) denotes the azimuthal angle 
with respect to the reaction plane. In this frame, w„ can be expressed as a function of the one-particle momentum 
distribution dNj/(fip for a particle of type j (in this paper, we consider pions and protons): 

cos ntp ^ dip 

VniPT,y,j) = (cOSn0} = ,27r ^ (1) 



where the brackets denote an average value over many events with approximately the same impact parameter, and 
Pt and y are the transverse momentum and rapidity of the particle. Since the system is symmetric with respect to 
the reaction plane for spherical nuclei, {smn(j)) vanishes. The goal of the flow analysis is to extract u„ from the data. 
The coefficients vi and V2 are usually called directed and elliptic flow, respectively p^ . 

Since the reaction plane is not known experimentally, one can only measure relative azimuthal angles between the 
outgoing particles. In particular, one measures the Fourier coefficients of the relative azimuthal distribution between 
two species of particles j and k 

dN,k , , , , 

cos 71(01 -92) 



C"'i(pTi,yi,j;PT2,2;2,fc) ^ (cosn(0i - ^2)) = '-^ w.r " • (2) 



d-'pid^P2 

dNjk 



d-'pid^P2 



The two-particle distribution can generally be expressed as the sum of an uncorrelated distribution and two-particle 
correlations: 

dNjk _ dN, dN, 

-75 - -J3 -Jt, (1 + <--jA;(Pl,P2j) (oj 

where Cjfe(pi, P2) is the two-particle connected correlation function. Using this decomposition, one can write c™°'*''"'°'^ 
as the sum of two terms: 

^measured {p^,, y, , j ; p^, , y,, k) ^ C«°- {PtI, yi , J PT2, y2, k) + c-"-«°- (p^l , ^1 , J ; PT2 , 2/2 , fc) (4) 

where the first term is due to flow: 

Cn°^{PTl,yi,j;PT2,y2,k) = Vn{pTl,yi,j)Vn{PT2,y2,k) (5) 

while the remaining term comes from direct two-particle correlations: 



C0S7l((/)i - 02)C-,fe(pi,P2' 



dNj dNk 



non-flow/ • i\ JJ d pi d P2 / r,\ 

C„ {pTuyi,j;PT2,y2,k) = jrj ^^^y- ^ ^ . (6) 



d^Pid'p2 

From the flow term (||), one can calculate the Fourier coefficient u„, up to a sign (see for details): 

, cfl°'-(pr,;/,j;P) 

Vn(PT,y, l) — =^ (7) 

v/c«°-(2?;2?) 

where V denotes the detector used to estimate the reaction plane, Cn(j>Ti,yi, j]T^) the average value of c„ over 
iPT2,y2,k) in 2?, and c„(23;2?) the average over both (pti,2/i, j) and {pT2,y2,k). 

Experimental analyses usually assume that the only azimuthal correlation between outgoing particles is due to their 
correlation with the reaction plane, i.e. they neglect c^on-Aow means that the published flow data are given by 
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Eq.(^ with cf°^ replaced by the total correlation c™'^'^*'^'''^'*; in turn, the total correlation can be reconstructed from 
the pubHshed flow data using Eq.(||). 

We shall see in section II] that c"°"^*°^ is typically of order where N denotes the number of outgoing particles. 
In a central Pb~Pb collision at SPS, N ^ 2500, so that one expects c^o°-fl°-«' ~ 4 x 10"''. This is certainly a weak 
effect. However, the total azimuthal correlation is also weak: since the directed and elliptic flow at SPS measured 
by NA49 are of the order of 3% §, one obtains from Eq.(||) (since it is assumed that all the measured correlation is 
from flow) 4°°^™'^°'^ ~ 9 X 10"'', only a factor of 2 larger than non-flow correlations. Therefore one cannot neglect the 
latter. In order to take them into account, one must estimate the various contributions to subtract them 

from c™^'*'^"'^'^^ in (^) to isolate c^°'^, and calculate the resulting Vn using (|^). In [Q, we applied this procedure to 
the NA49 data on pion flow taking into account correlations due to the HBT effect. In section [V, we consider in 
addition the effects of momentum conservation and resonance decays. 



III. SOURCES OF NON-FLOW CORRELATIONS 



Several physical effects yield two-particle correlations, which contribute to the non-flow correlation defined in (|^). 
In this section, we discuss momentum conservation, resonance decays, quantum HBT correlations, and final state 
interactions, and we estimate their respective contributions to c5^°'^~^°^. 



A. Transverse momentum conservation 



Global momentum conservation gives a back-to-back correlation between the produced particles. Since we are 
interested in azimuthal correlations, we consider only transverse momentum. Its contribution to Cjfe(pi,p2), which 
we denote by Cj^^^ (pi, P2), is calculated in Appendix |^ in the case where there is no other correlation between 
particles. The result is 



C^, (P1,P2)--^^^ (8) 



where the sum in the denominator runs over all the particles emitted in the collision. This correlation is clearly of 
order 1/N. The Fourier coefficient c^''^ is obtained by inserting (|^) in Eq.(j|). Neglecting azimuthal anisotropics of 
the one-particle distribution, which are of the order of a few percent , one finds a result which is non- vanishing only 
for n — 1: 

SpT / ■ i\ PT1PT2 

Ci (PTi,yi,j;pT2,2/2,fc) = --^^^- (9) 

Since c^^^ vanishes for n ^ 1, the effects of momentum conservation must be taken into account only in the mea- 
surement of directed flow. Eq.(|^) allows to calculate the correlation arising from momentum conservation, as soon 
as the value of the denominator (X]Pt) is known. This quantity is evaluated in Appendix ^ in the case of Pb-Pb 
collisions at SPS. The value (|B3|) is probably overestimated, which means that we underestimate the correction due 



to momentum conservation in the numerical calculations of section IV. Eq.(|9|) shows that the azimuthal correlations 
between two particles due to momentum conservation are independent of the rapidities yi and j/2 and of the particle 
species j and fc, but increase linearly with transverse momentum. The effect is therefore stronger for heavier particles, 
which have larger transverse momenta, as for instance protons. The centrality dependence is also easily determined: 
since transverse momentum spectra depend weakly on the centrality [p5 26|, (J^Pt) scales like the multiplicity N, so 



that cf^^ scales like 1/N. Note, finally, that momentum conservation does not contribute if the window used for the 
reaction plane determination is symmetric around mid-rapidity 



B. Resonance decays 



We now discuss the azimuthal correlation between the decay products of a resonance. It is known that many pions 
originate from p tttt decays, and that many nucleons which are excited into A resonances p7| , p8t , which then decay 
into a nuclcon and a pion. When pions are used to determine the reaction plane, p (resp. A) decays must be taken 
into account in the analysis of the pion (resp. proton) flow; when protons are used, A decays must be taken into 
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account in the analysis of the pion flow. Resonances with higher masses are less abundant and will not be considered 
in this paper; however, they may also give rise to sizeable correlations. 

The contribution of a resonance decay, say, A pn, to is the product of two factors: the first factor is the 

average value of cosn((/)i — 4)2) where 0i and (f)2 are the azimuthal angles of a proton and a pion originating from the 
same A; the second factor is the probability that a given proton and a g iven pion originate from the same A. This 
probability is hard to estimate. There are recent measurements of and A++ ||2^ multiplicities, but this is not 

sufRcient: indeed, due to the short resonance lifetimes, the decay products may interact again after they have been 
produced, in a way which depends on the detailed collision dynamics. Microscopic models of heavy ion collisions (see 
pO| for a review) could be used to calculate correlations from resonance decays, but the discussion in this paper will 
remain at a more qualitative level. 

Generally, in a two-particle decay, the azimuthal correlation depends on the relative magnitude of the velocity V 
of the decaying particle and the velocities Vi and V2 of the products in the rest frame of the decaying particle: in 
particular, the azimuthal angles of the decay products are equal (0i — (f>2) in the limit where Vi = V2 = 0, while they 
are back-to-back (0i = 02 +''') if = 0. This is illustrated in Figj^, which displays the average value of (cos(0i — (1)2)) 
and (cos2((/)i — 02)) for the decays A — > pn and p —> tttt, as functions of the transverse momentum of the decaying 
resonance. 




2000 




2000 



FIG. 1. Average value of cos(0i — ^2) (full curve) and cos 2(01 — ^2) (dashed curve), where 0i and 02 are the azimuthal 



angles of the decay products of A 
particle in MeV/c. 



pTT (left) or p ^ TTTT (right), as a function of the transverse momentum of the decaying 



Assuming that the transverse momentum distribution is exponential in the transverse mass tot — ^/Pr^^n^ , i-e. 
dN/(PpT oc exp(— tot/2^) (see Eq.(Bl)), one can average the correlations displayed in Fig. ^ over pT- The result is 
shown in Fig.^j, as a function of the inverse slope parameter T. 
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FIG. 2. Same as Fig.l, as a function of the inverse slope parameter T of the resonance distribution in MeV. Left: A 
Right: p —* TTTT. 
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Taking for A the same value of T as for protons, i.e. T ~ 300 MeV, one sees that the contribution to ci is large 
and positive, while the contribution to C2 is close to zero. In the case of p decays, assuming that T varies linearly 
with the mass one may choose T ~ 270 MeV, which results in both ci and C2 slightly negative. 

In order to obtain the contribution of A decays to c'^°'^~^°'" , one must multiply the correlation displayed in Fig.(|l|) 
by the probability pA that a given pion and a given proton originate from the same A, as explained above. Let a 
denote the fraction of nucleons are excited into A resonances. Then, the probability pA is a/A'^ where N^^ is the 
number of pions. From isospin symmetry, one expects this probability to be the same for neutral and charged pions. 
Since most particles are pions, TV^ ~ iV. If a is of order unity, we obtain finally c^~*^'^ ~ The centrality 

dependence of this correlation is also in if the relative abundance of A is independent of centrality. Similar 
arguments hold for the correlations due to p decays. 

Note that a way to eliminate correlations from A decays experimentally would be to measure protons and pions in 
widely separated rapidity windows. Indeed, the maximum rapidity difference between the decay products of a A is 
1.5. In the case of the p tttt decay, the maximum rapidity difference between the outgoing pions is 3.3, so that the 
corresponding correlation is harder to eliminate. 

C. Small angle correlations 

Finally, there are correlations between particles with low relative velocities: quantum HBT correlations between 
identical particles, Coulomb interactions between charged particles, and strong interactions. 

Azimuthal correlations between identical particles due to the HBT effect were studied in detail in [|l^ . The Fourier 
coefficient c™"'"(pi, P2) is of order unity if the relative momentum |pi — P2I is smaller than h/R, where R is the size of 
the particle source. It depends weakly on n. The order of magnitude of c™"^ integrated over phase space is therefore 
c™""" ~ (h/Rp)^, where p denotes a typical momentum. In the case of a high temperature massless pion gas, p is of 
order T/c, while the particle density N/R^ is of order (T/hc)^, so that c™""" is again of order 1/iV One would 
expect the centrality dependence of c™"^ to follow the same behavior. However, the measured HBT radii vary more 
slowly than N^^^ ||3^. In collisions with sulphur projectiles, a dependence in i? oc was found with a ~ 0.2 Q. If 
the same behavior holds for Pb-Pb collisions, the dependence of azimuthal HBT correlations with centrality will be 
in N~'^-^. Finally, c™'^{'D) is proportional to the one-particle momentum distribution, hence it is larger at low pT- 

More generally, all particles with low relative velocities undergo final state interactions, strong and/or electrostatic 
pgf . One defines the invariant momentum as Qinv = 2pvi-c\, where p, = mim2/{mi + m2) is the reduced mass 
and v-cci = Y^l — m^m|/ (Pi ■ P2Y is the relative velocity of the two particles with 4-momenta Pi and P2. Then 
the correlations due to strong interactions are sizeable only in a limited Qinv range, independent of the size of 
the interaction region. For instance, proton-proton correlations, which are dominated by strong interactions, are 
maximum for Qinv — 20 MeV/c, and negligible for Qinv ^ 40 MeV/c, for any projectile and target On the other 
hand, the range over which Coulomb correlations are significant depends on the system size. For Pb-Pb collisions at 
SPS, p — 7r+ correlations, which are dominated by the Coulomb repulsion, extend up to Qinv — 50 MeV/c Q. Both 
Coulomb and strong interactions give rise to small angle azimuthal correlations, similar to those due to the HBT 
effect. Thus, the resulting contribution to c"°"~''°™, denoted by cj^^\ is expected to be weakly dependent on the order 
n of the Fourier component, and c^^^{pT,y,'D) is roughly proportional to the one-particle distribution. 

The centrality dependence differs for Coulomb and strong interactions. For strong interactions, the strength of the 
correlation at a given Qinv varies like 1/V, where V is the volume of the system ||l^]. This is due to the fact that 
the range of the interaction is much smaller than the size of the system. Thus the centrality dependence follows a 
1/A^ behavior. For Coulomb interactions, which are long-ranged, the N dependence is weaker. We do not attempt 
to evaluate this dependence. By analogy with the Coulomb potential, it might vary like 1/L, with L the size of the 
interaction region, i.e. like 1/7V^/^. 

Both HBT and final state interactions can be eliminated by studying the correlation between particles widely 
separated in phase space, so that their invariant relative momentum is much larger than, typically, 50 MeV/c. 

D. Summary 

The properties of the three main types of two-particle correlations studied above are summarized in Table I. Note 
that while the average over all phase space of the three is of order 1/A^, they can be much larger in definite regions of 
phase space, for example at low relative momentum for HBT correlations. Note also that this catalogue of two-particle 
correlations may not be exhaustive. 
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As we have seen, two-particle correlations are generally of order 1 /iV, much smaller than unity. We therefore assume 
that the contributions of the various effects can simply be added. If one studies the correlation between identical 
pions, for instance, there are contributions to cj^on-flow fj-Qj^j momentum conservation, denoted by c^^^ , from p ^ tttt 
decays, denoted by c^^, from HBT correlations, denoted by c™"^, and from final state interactions, denoted by c^^^- 
The resulting non-flow correlation is 



^non— flow 



(10) 





Momentum conservation 


Resonance decays 


Small angle correlations 


Which pairs? 


any 


I2/1 -2/2I < 1.5 (A) 
|yi-y2| <3.3 (p) 


low relative velocities 
identical particles (HBT) 


n dependence 


n — 1 only 


non-trivial 


weak 


centrality dependence 


1/N 


1/N 


l/N"-*" ? (HBT) 
1/iV (strong) 
l/A^i/^ 7 (Coulomb) 


Pt range 


high Pt 


non-trivial 


low Pt 


under control? 


calculated 


simulations, microscopic models 


measured 



TABLE I: Properties of the various sources of two-particle azimuthal correlations 



IV. APPLICATION TO SPS DATA 



The most detailed flow analyses in Pb-Pb collisions at SPS were performed by the NA49 and WA98 collaborations. 
NA49 estimates the reaction plane using pions in the forward hemisphere, while WA98 uses protons and fragments 
close to the target rapidity. In this section, we discuss to what extent non-flow correlations may contribute to the 
measured azimuthal correlations, which are assumed to be solely due to flow in both experiments. In the case of the 
NA49 data, we shall explicitly subtract the various non-flow correlations, following the method outlined in section ||, 
and see what the "true" flow might be. Our discussion of WA98 data will remain at a semi-quantitative level due to 
the complex acceptance of the Plastic Ball detector used for the flow analysis. 



A. NA49 Data 



The NA49 experiment at CERN measures the directed and elliptic flow of pions and protons in Pb-Pb collisions at 
158 GeV per nucleon Charged pions are used to estimate the reaction plane, in the kinematic window 4 < y < 6 
and < PT < 600 MeV/c for directed flow, 3.5 < y < 5 and < pt < 2000 MeV/c for elhptic flow The azimuthal 
distribution of identified particles (protons and charged pions) is then measured with respect to this reaction plane. 
We now show the effect of taking into account the various sources of non-flow correlations discussed in the previous 
section. 

We begin with the directed flow of protons. It is obtained from from Eq.(^, where T> refers to charged pions in the 
above mentioned kinematic window. The numerator in this equation corresponds to the correlations between protons 
and the pions in V, while the denominator corresponds to the correlations between the pions in P. The arbitrary sign 
in Eq.(|7|) is chosen negative by NA49, in order to obtain a positive vi for protons (see Fig.||), as usually assumed at 
high energies in the forward rapidity region [ |l5| . 

The correction to the denominator in Eq.(M) from non-flow correlations will be discussed below, when we discuss 
the pion flow. It is a small correction, of a few percent. The numerator of Eq.(^ gets contributions from momentum 
conservation, A — > pn decays, and final state interactions. 

We first discuss momentum conservation. The corresponding azimuthal correlation is given by Eqs.(^) and ( |B3| ). 
Since it is negative, once it is subtracted from the measured correlation (^), this gives a positive correction to the 
correlations due to flow c^°". Because of the negative sign in (^, the correction to vi is finally negative. As can be 
seen in FigJ|, the correction is so large that the sign of the proton directed fiow is now negative for almost all pT- 
Now, the sign of the proton directed flow is particularly important since it is used to fix the arbitrary sign in Eq. (^ , 
which is chosen so that ui > in the forward rapidity region for protons. The surprising result is that the directed 
fiows of pions and protons now have the same sign. This is contrary to the theoretical expectation that nucleons and 
pions fiow in opposite directions, because of pion rescattering |]3^ , p9| , and to other measurements at AGS Q and at 
SPS 0. 
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Since the correlation between pions and protons from A — > p7r decays is positive (see Fig.Q), it may restore the 
original sign of the proton directed flow vi. As explained in the previous section, an accurate computation of this 
correlation would require a microscopic model. We present in Fig.^ a crude estimate of the effect, based on a simple 
Monte-Carlo calculation taking into account the acceptance windows of NA49 for protons and pions. We do not 
take into the possible rescatterings of the decay products. Our assumptions are the following: 50% of the outgoing 
protons originate from A decays, as suggested by simulations based on the UrQMD model [g8|. We take into account 
the fact that only 2/3 of A decays yield a charged pion. The mass distribution of A's follows a lorentzian, their 
rapidity distribution is flat between and 5.8 (i.e. between the target and projectile rapidities), and their transverse 
momentum distribution is exponential in tot (see Eq.(Bl)) with an inverse slope parameter T — 300 MeV, as for 
protons. One sees in Fig.§ that correlations from A decays can indeed restore the positive sign for the proton directed 
flow, since their contribution is almost as large as that from momentum conservation in absolute value. 
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FIG. 3. Directed flow vi of protons, integrated over the range 3 < y < 6, as a function of the transverse momentum pr, 
measured by NA49 (open squares), taking into account momentum conservation only (with error bars), and with a simulation 
of correlations from A decays (filled squares). Error bars are taken from experiment. 



In principle, one should also take into account final state interactions. Since the pions used to determine the reaction 
plane are approximately half 7r+ and half tt~ , Coulomb effects cancel. Correlations between pions and protons due 
to strong interactions are small and should not contribute significantly to c™''^''"'''''^. 

We now turn to the directed flow of charged pions. Since pions are also used to estimate the reaction plane, pion- 
pion correlations are involved in both the numerator and the denominator of Eq.(j^). Following Eq.([To|), we must 
take into account momentum conservation, p decays, the HBT effect, and final state interactions. We neglect final 
state interactions for the same reasons as for protons. According to Fig.||, the correlation from p decay is small. We 
also neglect it. Three sets of points are displayed in Fig.^ NA49 data without correction, with correction for HBT 
correlations (taken from p^]) and with correction for HBT and momentum conservation. While the correction due 
to the HBT effect is important at low p^, momentum conservation changes the pion flow at high pT', the latter effect 
is less important than for protons, because pions have lower transverse momenta. However, the positive pion flow 
observed by NA49 at high pT seems to be explained by momentum conservation. 

Let us say a few words of elliptic flow. Here, momentum conservation does not contribute, as discussed in the 
previ ous section. For protons, A decays probably give a negligible contribution to c™'^'**'"''^ as explained in section 
HIB. For pions, HBT correlations significantly redu ce elliptic flow at low pT, as shown in jlj. Here, p decays should 
be taken into account. As discussed in section III B| , they give a negative correlation, i.e. opposite to the correlation 
due to flow. Therefore, taking p decays into account would lead to a (slight) increase of the pion V2- 

Finally, the centrality dependence of the directed and elliptic flow of pions has been studied recently [Q. It is 
striking to note that both increase continuously as the reaction becomes more peripheral, which is the behavior 



7 



expected for non-flow correlations. We have shown that at least a fraction of the measured correlation, which is 
interpreted as flow, is due to non-flow correlations. Even if it is only a small fraction, it will increase with impact 
parameter and eventually dominate for the most peripheral collisions. It is thus likely that the measurements of flow 
at large impact parameter are contaminated by non-flow phenomena. However, none of the above mechanisms can 
explain the observed y dependence of directed flow, which increases (in absolute value) up to the projectile rapidity 
for peripheral collisions. 
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FIG. 4. Directed flow ui of pions, integrated between 4 < y < 5, as a function of the transverse momentum pt in MeV/c, 
measured by NA49 (open squares), after subtraction of HBT correlations (filled squares), and after subtraction of both HBT 
and transverse momentum conservation correlations (with error bars). 

B. WA98 Data 

We shall now discuss the results obtained by the WA98 experiment at CERN for the directed flow of protons and 
positive pions There, the reaction plane is estimated using protons and fragments detected by the Plastic Ball, 
in the target rapidity region —0.6 < y < 0.3. One then measures the azimuthal angle of a positive pion or a proton 
detected in the Plastic Ball with respect to this estimated plane. If the particle is a proton, it is not used in the 
determination of the reaction plane to avoid autocorrelations. Two points closer to midrapidity are obtained using 
the tracking system, for which the flow is zero within error bars; we shall not consider these two points here. Contrary 
to the NA49 analysis where pions were used to estimate the reaction plane, the flow of protons (resp. pions) is now 
determined from proton-proton (resp. pion-proton) azimuthal correlations. The value of vi is still obtained from 
Eq.(^, where the arbitary sign is again chosen negative, so that the directed flow of protons be negative, as it should 
in the backward rapidity region. On the other hand, the directed flow of pions is positive, i.e. pions flow opposite to 
the protons. 

WA98 measures the directed flow vi in the target fragmentation region. Since vi vanishes at midrapidity, it is 
expected to be larger there than in the more central rapidity region, where NA49 is working. Indeed, the measured 
values of vi are larger for WA98 than for NA49. For semi-central collisions, |wi| is about 0.1 for both pions and 
protons. Following Eq.(||), this corresponds to an azimuthal correlation | | ~ 10~^, either between pions and 
protons or between protons, an order of magnitude larger than for NA49. 

The values of vi are measured as a function of centrality. For protons, —vi reaches a maximum of 0.12 for an 
impact parameter 6 ~ 8 fm, while for pions vi increases with b up to very peripheral collisions, where it reaches values 
as high as 0.3. The latter centrality dependence is reminiscent of the 1/N behavior of non-flow correlations, which 
suggests that part of the measured flow is spurious. We thus review the various non-flow correlations listed in section 
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Ill to evaluate their contributions. We shall see that final state interactions are the only significant ones, and that 
they may explain the pion flow at large impact parameter. 

Momentum conservation is less important than for the NA49 experiment. The reason is that the particles detected 
in the Plastic Ball have low transverse momenta For protons, px is typically 400 MeV, resulting in cf^^ of the 
order of —4 x 10^^ for the proton-proton correlation, following Eq.(|^). This is much smaller than q^^^^^^^'^ ^ 10^^. 
For pions, transverse momenta are even smaller. In both cases, the contribution of transverse momentum conservation 
to c5'™-f>°"' is negligible. 

Resonance decays will contribute to the pion flow. The sign of the correlation due to A — > pir depends on the inverse 
slope parameter T, as can be seen in Fig.^. This quantity is smaller in the fragmentation region than in the central 
rapidity region, but probably still larger than 100 MeV. Then, the corresponding pion-proton correlation should be 
positive. It is thus opposite to the correlation measured by WA98, and cannot explain the observed pion flow. 

Finally, let us discuss the small angle correlations due to the HBT effect and final state interactions. We shall 
first summarize the results obtained by the NA49 experiment for two-particle small angle correlations in the central 
rapidity region. Then we extrapolate these results to the fragmentation region covered by the WA98 Plastic Ball in 
order to discuss the correlations measured by WA98. 

The NA49 collaboration has measured both proton-proton and proton-7r+ |36j small angle correlations for 
central collisions. Two-proton correlations, resulting from the interplay of the HBT effect and strong and Coulomb 
interactions, were found to be positive, which means that they are dominated by the attractive strong interaction. 
They peak at a relative momentum between particles of 20 MeV/c, where they reach 25%, and are no longer significant 
as soon as |pi — P2I > 40 MeV/c. Small angle correlations between positive pions and protons are negative, and 
typically of the order of —10% for invariant momenta smaller than 40 MeV/c, i.e. Iptit — (m-ir/m-p)PTp| < 40 MeV/c 
and y-K — Vp < AQ/ttIt^ ~ 0.3. 

Extrapolating these results to the fragmentation region, we can estimate the azimuthal correlation cf^^ for WA98. 
In the case of proton-proton correlations, we need to evaluate the fraction of proton pairs having relative momenta less 
than 40 MeV/c, i.e. with |pti - PtjK = 40 MeV/c and 2/1 - 2/2 < Ap/mr ^ 0.04, with mr ^ 1000 MeV. Since 
the acceptance window is typically Q pf?" = 200 < pr < Pt""" = 600 MeV/c and y™'" = -0.4 <y < y""''' = 0.1, 
this fraction 



-2j< 
11 

4f Ap3 



5 X 10"". (11) 



With a correlation strength of 20%, this yields a value of cj^^ of the order of 10"*, much smaller than the measured 
correlation 2™°''''"''°'^. Thus, the measurement of the proton flow is not influenced by small angle correlations. 

In the case of pion-proton correlations, we flrst note that the velocity windows of pions and protons of the WA98 



Plastic Ball overlap |41|, so that small angle correlations between pions and protons should be taken into account. In 
fact, the rapidity range is almost the same for pions and protons, while the transverse momentum range is roughly 
scaled by a factor mj^/mp, so that the acceptance windows roughly coincide in velocity space. Then, for a given pion, 
the probability that a proton seen in the Plastic Ball has an invariant relative momentum Ipttt — i'mTv / 'm-p)PTp\ < 
Ap = 40 MeV/c is given by ( pl| ) with Ap replaced by {mp/m^)Ap. This gives a factor {mp/mT^Y — 200, compared 
to Eq. (|Tl|) , so that the fraction is of the order of 10~^. Using the correlation strength measured by NA49, this 
gives a correlation cf^'(2?) of order — 10~^, comparable to the measured correlation. On the one hand, this value is 
overestimated since the acceptance windows for pions and protons do not exactly coincide in velocity space. On the 
other hand, it is underestimated because the correlations are measured by NA49 for central collisions, and they are 
even stronger for peripheral collisions. Therefore it is likely that a large fraction of the 7r+ directed flow seen by WA98 
is due to the repulsive Coulomb interaction between positive pions and protons. 



V. DISCUSSION 



We have shown that the assumption that all azimuthal correlations between particles are due to flow is no longer 
valid at SPS energies, where other sources of correlations become of the same order of magnitude due to the smallness 
of the flow. All the effects studied above, momentum conservation, resonance decays, HBT correlations, final state 
interactions, turn out to give important correlations. In particular, we have shown that the effect of momentum 
conservation is so large that it changes the sign of the proton directed flow measured by NA49: this sign is particularly 
important since it is used to deflne the orientation of the reaction plane. 

The methods currently used to analyse the flow can be modified in order to take into account these additional 
correlations and to subtract them from the measured correlations, as explained in section ^ Correlations from 
momentum conservation can be calculated. They are governed by only one unknown quantity, the sum over all 
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particles of squared transverse momenta (X^-Pt)- Correlations from the decays of short-lived resonances cannot be 
evaluated so easily. In this paper, we have only given semi-quantitative estimates. More accurate calculations would 
require a microscopic model of the collision. Finally, correlations which affect particles with low relative velocities, i.e. 
those due to the HBT effect and to final state (Coulomb, strong) interactions can either be measured independently, 
or eliminated by measuring correlations between particles in widely separated rapidity and/or transverse velocity 
windows. 

However, our list of correlations may not be exhaustive, and one cannot exclude that other sources exist which 
are of the same order of magnitude as those studied here. There is no systematic way to separate the flow and the 
non-flow contribution to the measured azimuthal correlation (^). One could try to use the factorization property of 
Eq.(||) to isolate the correlation due to flow. However, non-flow correlations may also factorize, as is the case for the 
correlation due to momentum conservation (|^). 

A way to circumvent these difficulties would be to use the remarkable property of two-particle correlations, that 
they vary with the multiplicity N like 1/N (with the possible exception of Coulomb interactions and, to a lesser 
extent, HBT correlations). This fact has been noted in section [II for the various correlations studied, but holds 
more generally for any extensive system, and would probably still be true for other types of correlations than those 



listed in section HI. This 1/A^ behavior determines the centrality dependence of non-flow two-particle correlations. 
Since the flow is known to vanish for central collisions by symmetry, one could measure (and then subtract) non-flow 
correlations by measuring them for central collisions and then assuming that they vary like 1/N. Unfortunately, it is 
impossible to select "true" central collisions (with vanishing impact parameter) experimentally, so that this procedure 
may not be accurate. 

Future heavy ion experiments at RHIC and LHC are likely to be affected in the same way by non-flow correlations. 
At these energies, directed flow is expected to be weaker than at SPS, maybe undetectable. On the other hand, 
elliptic flow should be at least of the same magnitude as at SPS. Pions should be used to determine the reaction 
plane since they are by far the most abundant. As the multiplicity will be higher than at SPS, non-flow correlations 
will be smaller, since they are proportional to 1/N. However, since the detectors will cover only a limited interval in 
rapidity, effects of p decays, HBT correlations and final state interactions may be important and should be taken into 
account in the measurement of the pion elliptic flow. Hard processes could give an additional contribution to non- 
flow correlations: there is a strong azimuthal correlation between jet fragments Clearly, more work is needed, 
both theoretical and experimental, in order to obtain accurate measurements of flow observables at ultrarelativistic 
energies. Alternative methods, based on multiparticle correlations, are currently under study [H. 
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APPENDIX A: CORRELATIONS FROM MOMENTUM CONSERVATION 



We denote by pi • • • pat the momenta of the N particles emitted in a heavy ion coUision and by pxi ■ • • Ptat their 
transverse components. Since we are interested in azimuthal correlations, we consider only transverse momentum 
conservation: pxi -I- ... -I- Pttv = 0- We assume for simplicity that this is the only correlation between the momenta 
of the outgoing particles. We want to calculate the distribution of pi, . . . , p^, with k < N, defined as: 

fciPi,...,Pk)^^ JT^^ (Al) 

J S^PTi + ■■■+ PTn) n {fiP^)d'P^) 

2—1 
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where /(p) denotes the single particle normalized transverse momentum distribution. The average transverse mo- 
mentum is naturally assumed to be zero: 



(PT> = J PT/(p)rf'p = 0. 



(A2) 



In order to calculate the distribution fc defined in Eq.(Al), we make use of the following lemma: the sum of 
M uncorrelated momenta Pt = J2i=i PTi ^ gaussian distribution if M is large, according to the central limit 
theorem: 



» / M \ M 

Fm(Pt) ^ / -Pt + 5]pt. n (/(pO^'pO 

\ i=i J i=i 

1 / Pt^\ 

= — exp — . 



(A3) 



In writing this equation, we have assumed that the transverse momentum distribution is isotropic in the transverse 
plane, i.e. we have neglected the azimuthal asymmetries (|l]), which are of the order of a few percent. The width of 
the gaussian is 



m{pI). 



Eqs.(|A3|) and (|AJ) allow to write (|AlD as 



(A4) 



/c(Pl 



■,pfe)= (^n/(p.)^ 



Fw(0) 



N 



N -k 



exp 



(Ept.)' 



V 



For k = 1, Eq.(A5) gives the corrected one-particle distribution 



/c(p) = /(P) 



1 ^ Pt' a 



where we have expanded the correction to leading order in Similarly, one obtains for k—2 

2 (PT1+PT2! 



/c(Pl,P2)=/(Pl)/(P2) 1 + 



N N{p\) 

The contribution of momentum conservation to the correlation function (0) is thus 



r^PT(^ „ A /c(Pl.P2) -1 2pTi ■ PT2 

/c(Pl)/c(p2) iV(PT) 



(A5) 



(A6) 



(A7) 



(A8) 



Similar results were obtained in llq]. Appendix B. 



APPENDIX B: ESTIMATE OF ^^(^t) AT SPS 

In a central Pb-Pb collision at SPS, about 680 negatively charged particles arc emitted |45j. If we assume that 
these are mostly tt^ and with a ratio {K~) / (it) ~ 0.09, where (tt) = ((7r+) -I- (7r^))/2 |2^], there are some 625 tt^ 
and 55 K~ . Taking into account the various kaon-to-7r ratios, as well as the expected isospin symmetry between pion 
species, this gives a total of 280 kaons and 1875 pions. With these particles come the 416 nucleons (we assume there 
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is no fragment), giving about 2500 particles in all. We shall neglect the contributions to YliiPr) from other particles 
as e.g. antiprotons. 

In order to calculate the average value of the squared transverse momentum {p^) for a particular particle species, 
one needs its distribution in transverse momentum and rapidity. The px- and y-dependences may factorize, leaving 
a normalized pT distribution which is parametrized by 



dN e"/^ 

■ exp 



(P^T 27rT(TO + T) 



(-^)^ (Bl) 



where rriT = \/ vn? + p\. is the transverse mass, while the inverse slope parameter T depends on the particle species 
and the energy of the collision. This parametrization can be used for the pions, nucleons and kaons at SPS, with 
= 180 MeV, Ttv = 300 MeV and Tk — 215 MeV respectively [Q. These numerical values fit transverse momentum 
spectra in the central rapidity region. Away from mid-rapidity, transverse momenta are somewhat smaller, so that 



we overestimate Ep|,. The spectrum (Bl) yields an average (p^) given by 



1 + 3— + 3— 

(p?,) = 2mT " y (B2) 

Substituting in Eq. (|B2| ) the values of the inverse slope parameter T given above, and multiplying the resulting 
average squared transverse momentum by the multiplicity of the particle species, one obtains for a central collision 
at SPS. 



E 



- 930 GeV^ (B3) 



As mentioned above, this is somewhat overestimated since we have neglected the decrease of inverse slopes away from 
mid-rapidity. 

For a non-central collision, we assume that the sum of p^ scales with the multiplicity, i.e. we neglect the centrality 
dependence of inverse slopes and particle ratios. Therefore, in the case of semi-central collisions used by NA49 for 



flow studies, with only 40 to 55% of the maximum multiplicity M, the value (B3) must be divided by 2 
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